Spin density wave formation in graphene 
facilitated by the in-plane magnetic field 
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We suggest that by applying a magnetic field lying in the plane of graphene layer one may facilitate 
an excitonic condensation of electron-hole pairs with opposite spins and chiralities. The provided 
calculations yield a conservative estimate for the transition temperature Tc 0.1 B. 



Charge Density Wave (CDW) and Spin Density Wave 
(SDW) formation in graphene is an intriguing, possibility 
which has been discussed in the literature Even 
more so since the reduced dimensionality of graphene will 
strongly affect the character of the corresponding ordered 
state probably making it critical. It was also pointed out 
that magnetic field perpendicular to the graphene lay- 
ers may facilitate formation of the CDW The 
CDW order parameter discussed in these papers estab- 
lishes different population densities on the two sublattices 
of graphene (we will call them u and v ones). This corre- 
sponds to a site-centered CDW. In this paper we suggest 
that by applying a magnetic field in the graphene plane 
or a Weiss exchange field 0] one can facilitate a formation 
of a bond-centered spin density wave. The corresponding 
order parameter has U(l)xZ2 symmetry and combines 
fermion operators with opposite spins and chiralities. 

The tight binding Hamiltonian is 



Hn = -t 



:(r)w(r + hi) + H.c. 
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where r belong to a triangular lattice (see Fig. 1). 
• V 





discuss this small deformation of the spectrum since it 
does not destroy the nesting between the bands of dif- 
ferent chirality and therefore will not have any adverse 
effect on the phenomena discussed in this paper. 

The Hamiltonian describing states close to the tips of 
the Dirac cones is 
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Here u, v operators annihilate electrons on different sub- 
lattices and k = kx + iky, Q = (1, l/\/3)27r/a\/3. Since 
magnetic field B lies in the plane, it affects only spin (the 
same effect can probably be achieved by the exchange 
field when one brings graphene sample in a close proxim- 
ity to a ferromagnet). The in-plane magnetic field leads 
to Zeeman splitting of the bands (see Fig. 2). The pre- 





FIG. 1: Crystalline lattice of graphene. 



FIG. 2: Graphene bands split by the in-plane magnetic field. 



Tunneling on next nearest neighbours generates the 
so-called trigonal warping of the spectrum. We will not dominant interaction is the Coulomb one: 
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y=i J2 [u+{k' + q + TQ)u^{k + q + TQ)+v+{k' + q + TQ)v^ik + q + TQ)]-'^''''' 



k.k' ,q:T — :tl 

X (k + t'Q)u^> [k' + t'Q) + v+ {k + T'Q)va, [k' + t'Q)] 

I 



|k-k' 



(3) 



The bare Coulomb interaction in graphene is fairly strong 
/hv but it is screened: 



V{q) 



\q\ - 27re2iVn(w, g) 
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where A'^ = 4 is the number of fermion species (valley 
and spin). Assuming that the self-consistent interaction 
is weak we can get an estimate of this interaction by 
substituting there the polarization loop Ti{Lj^q) for bare 
electrons. Then the effective interaction can be approxi- 
mated as 



V{q)^l^WF)]-' 
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where pi^F) — kp/2-K and C{N) ^ 1 is constant (accord- 
ing to 8 C(oo) =7r/4). 

The energetics selects order parameters where elec- 
trons pair with holes of the opposite spin and chirality. 
There are two complex order parameters: 

^l,{k) = Y.^u,^^' + - k'l)t'^.(k' - Q)), 

k' 

A^3(fc) = E("--(k' + Q)^(|k - k'|)«+(k' - Q))(6) 



The self-consistency conditions are 



Ai4(fc) 



(9) 

/1 2 
J^V{\k - p|) [nn,44(p)Ai4(p) + ni2,34(p)A23(p)] 

A23(fc) = (10) 

d p 



where 



rF(|k - p|) [n22,33(p)Ai4(p) -I- n2l,43(p)A23(p)] 



nn,44(p) = n22,33(p) = TY,G\l{^n.p)Gi^{u:n.p) 
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ni2.34(p) = n21,43(p) = T^G\2{uJn,p)G^4{uJn,p) 



As we said, these equations yield Ai4(fc) = A23(fc) — 
A(fc). Summing over the Matsubara frequencies and in- 
tegrating over angles we get the following equation: 



As we shall see, on the mean field level these order 
parameters have equal amplitudes independent of spin: 
A^4(fc) = Af3(fc) = A{k). Plugging §ii into ^ we ob- 
tain the mean field Hamiltonian: 
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The spectrum is 

El = |A(fc)|2+t;2(|fc|±B)2 



(8) 
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A(fc) 



In integrating over angles we will adopt the same approx- 
imation as in Q and will replace the resulting elliptic 
function K{k/p) by 1. In the limit of iV — > oo these 
equations yield Tc ^ Bexp{—N). We believe that us- 
ing this formula for A'^ = 4 substantially underestimates 
Tc. The approximate equations for the mean field critical 
temperature P* — B/2Tc » 1 is (here x =p/B): 



A(^) = T^r 



r tanh[(y + 1)(3*] tanh[(y - 1)13* 
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We found a good interpolation formula will split the in- 



terval of integration into (0,A) and {A, oo) where ^ > 1. 
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At X < A the interaction will be approximated as con- 
stant, and at a; > ^ as {Cx)~^. We look for the solution 
of this equation as 



A(x) 



(2; + 1)^ 



where 7 = 1/2 + ^1/4- l/CN and 

T,«Bexp[-(iV-l/7)] 



(13) 



(14) 



Assuming that C(4) « 1 we get 7(4) w 1/2 and the 
estimate ~ lO-^S. 

Few words about the effects of reduced dimension- 
ality. As is well known, in two dimensions continuos 
U(l) symmetry cannot be spontaneously broken. In this 
case the second order phase transition is replaced by the 
Berezinskii-Kosterlitz-Thouless (BKT) one and the low 
temperature phase is critical. The BKT termperature is 
determined by the equation 



BKT 



BKT) 



(15) 



where p(T) is the temperature dependent stiffness. p(T) 
vanishes at Tc and reaches some constant value at T = 0. 
Since the system is not strongly interacting, p(0) >> Tc- 
Therefore the BKT transition temperature occurs close 



to the mean field and one should expect to see the effects 
of fluctuations only in a fairly narrow region around the 
transition temperature. 

Our estimate of the transition temperature indicates 
that it is not too low and for realistic magnetic fields 
S ^ 10 — 20 T probably lies in the interval 1 K or even 
larger. In a hypothetical construction a graphene layer is 
brought into a microscopic contact with a ferromagnetic 
film; in this case the Weiss fields can probably be made 
higher by an order of magnitude. 
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